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Abstract 

Character expansion developed in real space renormalization group (RSRG) approach 
is applied to U(l) lattice gauge theory with #-term in 2 dimensions. Topological charge 
distribution P(Q) is shown to be of Gaussian form at any /^(inverse coupling constant). 
The partition function Z{9) at large volume is shown to be given by the elliptic theta 
function. It provides the information of the zeros of partition function as an analytic 
function of £ = e 16 (6 = theta parameter). These partition function zeros lead to the 
phase transition at 9 = n. Analytical results will be compared with the MC simulation 
results. In MC simulation, we adopt (i) "set method" and (ii) "trial function method" . 
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1. Introduction 



The numerical approach to the system with (9-term suffers from the difficulty, because 
the Boltzmann weight is complex valued and cannot be adopted as the probability weight 
directly. The character expansion method developed in the real space renormalization 
group approach is not affected by this difficulty @ 1|. It plays quite powerful role to 
understand the gauge systems with #-term 1| f|] • The calculation is based on the 
orthonormal property of the irreducible group characters. It can be performed in the 
system with #-term without any difficulty 

Moreover the existence of the #-term leads to some new features || [TO | [IT ||12|| [IB 



14] [15] in the nature of partition function and leads to new physical situation ||1 2\\ which 



are absent in the system without 0-term [|l6|| . 



We present in this paper the character expansion for the U(l) gauge theory with 
(9-term in two dimensions. The action is assumed to be the sum of Wilson real action 
with real (inverse) coupling (3 and standard imaginary #-term action with the coupling 
a = 9/2n. The expression for the partition function at any /3, i.e., in both strong and 
weak coupling regions, will be provided and it leads to an expression for topological charge 
(Q) distribution (P(Q)) at any (3. 

We obtain Gaussian distribution, P(Q) oc exp(— Ky (/3)Q 2 ), where Ky((3) is a constant 
determined by the value of real coupling constant (3 and the size of the volume V. The 
concrete expression of Ky(f3) will be given. In the weak coupling regions, it can be expressed 
as 

k v ((3)=C([3)/V<x(3/V, (1.1) 
and C((3) approaches x (3 in (3 — > oo limit ( ~ 19.7). In the strong coupling limit, 

C(0) = 6. (1.2) 

This result will be compared with numerical analysis. 

Here we note about the numerical approach[|T7| . As mentioned above, the Boltzmann 
factor for the action composed of both real and imaginary parts can not be adopted as 
probability weight. So we use the probability weight given by the real action only and 
calculate the topological charge distribution P(Q). Then the Fourier series of P{Q)e l6< ^ 
gives us the partition function. 

The comparison of P(Q) given by the character expansion method with that by nu- 
merical simulation will be given. We see quite good agreement between these approaches 
in all regions of the real coupling (3. 

In the MC approach to obtain P(Q), there is a technical difficulty related to the Q 
dependence of P(Q). The function P(Q) is quite rapidly decreasing function of Q at large 
\Q\ and the numerical simulation of P(Q) in large range of \Q\, at once is impossible. The 
ways to escape this difficulty are [pT7| p8[1 , 
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(i) The MC simulation is performed in the sets of smaller \Q\ range. This is said "set 
method " |T!| and afterwards the results in these sets are gathered by adjusting neigh- 
boring results. 

(ii) Even if the range is decomposed into sets, the obtained weight P(Q) still changes too 
rapidly in each set. So we use improved P(Q). Improved P(Q) is given by normalizing 
P{Q) by appropriate trial function of P(Q) which is obtained by foregoing small 
numerical simulation or by some theoretical prediction if we have it. 

In all regions of /3, the topological charge distribution is found to be expressed by 
Gaussian form P(Q) oc exp(— Ky ((3)Q 2 ) quite well and we have an interesting analytical 
expression for the partition function Zy(9). It is expressed as the third elliptic function 
$3 in the large volume limit (y — > oo). It leads to infinite series of zeros of Zy{&) on real 
negative axis of £(= e l9 ) plane, since $3 has the form of infinite product of linear function 
of C- The closest zero to the physical point 9 = 7r(£ = —1) suggests there is a first order 
phase transition [^TJ at 9 = 7r. The distance, Im6> = Ky(f3), of this closest zero to 
C = — 1 is proportional to 1/V at any (3 *' and an analytical expression for this distance 
is given as a function of j3. 



' Wiese obtained analytical estimate at (3 = 0. 
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2. Character Expansion 

(2.1) Character Expansion 

We will present the character expansion for the partition function of the U(l) gauge 
theory with #-term. The partition function is 

Z V {6) = f [dx l ]e^- s ^ +iv ^\ 



where x\ and x p (p = 1, V) denotes the link and the plaquette variable, respectively. 
The action is given by the Wilson real action and standard imaginary action, 

— s(x) + iv(x) = (3cosx + iax (2-1) 

where x = F i(F i means the field strength), /3 = 1/g 2 and a = 6/2n. Single plaquette 
contribution is 

F p (x) = exp{-s(L, x) + iv(L,x)} (2.2) 
The partition function for volume V is 

[d Xl ]F{x u ...,x v ) (2.3) 

-7T 

where F-function T is given by 

oo 

J 7 (x 1 ,...,x v ) = F(x 1 ,...,x v )8(x 1 + ... + x v -2ttQ). (2.4) 

Q= — oo 

In (2.4), Q denotes the topological charge (integer). By use of Poisson sum formula, T 
can be written as 

oo 

T{x u ...,x v )= J2 F(x 1 ,...,x v )e^ Xl+ - +x ^ (2.5) 

l= — oo 

The functional F(x\, ...,xy) is a product of plaquette contributions, 

F(x 1 ,...,x v ) = F p (x 1 )F p (x 2 ) F p (x v ) (2.6) 

The character expansion for F p is 



F p( x )= Yl F m(P,u)Xm(x), 

m= — oo 

Xm(x) = exp(zmx), m = integer, 



(2.7) 
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where Xm( x ) denotes the irreducible character specified by integer m. Then 

oo 

F(x u ...,x v )= F p ( Xl )F p (x 2 ) F p (x v )e u ^+-+^ 



l= — oo 
oo V 

Y X\^Frn k {P,Oi)x mk [Xk)Xl{Xk) 



l= — oo k=l m k 
oo V 



= Y II ^2 F m k W^ a )Xm k +l(Xk) 
l= — oo k=l m k 
oo V 

= Y IIE^-^'^mtW 

l= — oo k=l m fc 

The character expansion coefficient is0 

F m (/3,a)=^F 71 s (/3«_ n (a), 

n 

where character expansion of exp(— s) and exp(iv)(= exp(-iax)) are respectively, 

F n s (/?)=I n (/3), 

~ shi7r(a + /) 

7r(a + I) 

where I n ((3) denotes the modified Bessel function. It satisfies the relation 

F m+ i(/3, a) =F m (P,a + l) 

Then 

oo V 

F(xi,...,X V )= II Y Fm ^' a ~ l )Xm k {x k ) 

l= — oo fc=l m k 

Now we integrate all the inner links except for the outermost link variables. 

\\ dxi = i xmex \\ x ^J : {xi, xv) = 22 2jF m (/3, a - I) Xm{xi + — + %v) 

_7r l l m 

y = X\ + ... + xy = outermost link variable. 

where the orthogonality 

T dv 

(Xl + v)Xm 2 [X2 - v) = 5 miTn2 X mi {x\ + X 2 ) 

it 27T 



is repeatedly used. The partition function is 



Zv{9)= ^-J2J2 P rn(P,a-l) V Xm (y) 
J -« /7F I m 



= E J2( P rn(P, a - l)) v 5 m , (2.15) 

I m 

i 

where 



(2.2) Topological charge distribution 

Topological charge distribution appears in 



Z V (9)= P(Q)e l9Q , (2.17) 

Q = — oo 

and P(Q) is 

/•1/2 

P(Q) = / daZ v (e)e- l2naQ . (2.18) 

i-l/2 

Eqs.(2.15) and (2.17) lead to, 

•1/2 



P(Q) 



da V (F (/3,«-/)) y e- l27raQ 

-1/2 I= ~ 



/oo 
da(F ((3, a)) v e 
-oo 



(2.19) 



V -i27raQ 



where Fo(/?, a) is given by eq.(2.16) setting / = 0. This is an exact expression in f3. 

We will evaluate (2.19) in large V limit. The character coefficient Fo(/3,a) is an even 
function of a and peaked at a = 0. In large V limit, (Fq((3, a)) v is a function sharply 
peaked at a = 0. So we can evaluate its contribution by the Taylor expansion around 
a = 0. 

sin n(a + n) 



F (f3,a)= ^(PY- 



Ti(a + n) 



n= — oo 

oo 



T , n . simra n v-^ T . n . sin7r(a + n) 

W) + 2 )^ I n (J3) f, / 

7TQ; 7r(a + n) 



n=l 

00 



lo(0)(l - -a 2 ) - 2a 2 V InW-^- + °(« 4 ) 

n=l 



(2.20) 



7T 2 o 00 /_ i\n 

= /o(/3)exp(-G(/3)a 2 ), 
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with 



for any (3. 

Now we have 



cw = y + w)^ i " w -^-' (2 ' 21) 



P(Q) = 2 da exp(-G(f3)Va 2 )cos(2iraQ) 



o 



7T , -TV 2 2 



exp( ——— Q 2 ) (2.22) 



GV ^G(/3)y 



7T 



The quantity is defined by 



exp(-M/?)Q 2 ). 



7T 2 



«vW = (2-23) 



We define 



7T 2 



C(/?) = « y (/?)F= — (2.24) 

Topological charge distribution has the Gaussian form. The coefficient nv(P) has the 
"1/V" form for all (3. The coefficient C{(3) is plotted as a function of (3 in Fig.l. 



Fig- 1 

It is given, in the strong coupling limit, by 

G(0) = TP (2.25) 
C(0) = 6. 

In weak coupling regions (/3 > 1), it seems proportional to (3. The value of C((3)/{3, 
however, changes slightly with (3 (Fig. 2). 



Fig. 2 



(2.3) Partition function zeros and phase transition 

When P(Q) is given by Gaussian form (2.22), the partition function 

OO 

Z v {6) = Y,P{Q)^ iaQ = £ P(Q)C Q , (2.26) 

Q Q= — oo 
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has a simple factor ized form *) 



Zv(Q) = y^E ex P(-^(/3)Q 2 )e 27riaQ , 

I 00 2n— 1 

=v / ^n( i +^- i o( i +V )i 

n=l 



(2.27) 



where ( = e lfc/ = e 27 " a and 



?o=n( i -^ 2n )' 



n=l 



(2.28) 



q = exp(-Ky(/3)), (< 1). 



It can be written as 



where $3(1/, r) is the third elliptic theta function with q = exp(z7rr) and a = v. 

Important feature of (2.27) is that the partition function has infinite zeros (at V — > 00) 
on the real negative ( axis in complex ( plane. It shows that the system undergoes 
transition at 9 = tt. When V is finite, the closest zero of Zy{9) to ( = — 1(9 = n) is 
located at the distance Im# = kv(P) oc 1/V from £ = — 1 point. In V — > 00 limit, it 
approaches the ( = — 1 point and infinitely many zeros accumulate on the real negative 
axis of complex £ plane. 

Poisson sum formula applied to (2.17) gives 

exp(— KyQ 2 + 27riaQ) 

Q 

/oo 
d£ exp(-K V £ 2 + 2mai + 2%i£m) (2.30) 

' m 

Note that this is a periodic Gaussian form and that it shows clearly the phase structure 
as follows; For mo — 1/2 < a < mo + 1/2, the sum is dominated by m = mo and for 
mo + 1/2 < a < rriQ + 3/2 the sum is dominated by m = hiq + 1. It shows that there 
are infinitely many phase transitions at a = mo + 1/2, i.e., at 9 = it + 27rmo(mo = any 
integer). 

Eq. (2.30) can also be written as, 



*) The authors are grateful to Prof. Y. Yamada for valuable discussion on this point. 
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exp -— - )$>xp(-— m 2 -——m). (2.31) 

m 

The left hand side of (2.30) is -# 3 (z/, r). Eq.(2.31) is just the same form as the left 
hand side of (2.30) if we replace Ky ~" > K y = n2 / K v, a — » a' = ian/ny. Then we have 

M", r) = /I exp( — )0 3 (i/, r' ) (2.32) 

\ IT T 



where 



v = a, 
mr = —Ky 



(2.33) 

T = -1/T, 



z/ = u/t. 

Eq (2.33) is a kind of modular transformation. 



3. Comparison with numerical method 

When we adopt the action with 0-term, the Boltzmann factor 

exp(— s(x) + iv(x)) 

is a complex number and cannot be directly used as a probability weight. The partition 
function can be written 

Zv{0) = / [dxi] exp(- s(x p ) + ia x p ) 
J p p 

= E / [^] (Q) ex p(- E s M)^ iaQ (3.i) 



E p (^) 



where [ctez]^ is the integration measure restricted to topological charge Q. We first 
calculate the topological charge distribution P(Q) with real positive probability weight 
ex P(~ Yl s ( x p)) an d then we obtain Zy(0) through Fourier series. In order to perform 
numerical simulations to obtain P{Q), a technical problem appears because P(Q) is usually 
a quite rapidly decreasing function of Q and it drops by many decades at large Q. Two 
technical methods to avoid this difficulty will be used. 
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(i) "Set method": Wide range of Q is divided into many sets. Each set composed of 
[Q, Q + 1, Q + s — 1]. In the actual calculation we choose s = 4. Neighboring sets 
are adjusted at the overlapping Q so as to give P(Q) over these neighboring sets with 
the common normalization. Consider two neighboring sets [Q — (s — 1), • • • , Q] and 
[Q, Q + 1, Q + s — 1]. We equate the value of P(Q) at largest Q{= Q) of the lower 
set with that at the smallest Q(— Q) of the upper set. 

(ii) "Trial function method (Improved P(Q))": Even we divide into sets, the change of 
P(Q) in each set is still large. We define improved P(Q) dividing P(Q) by trial P(Q), 

" P t {QY 

P(Q) = P(Q)/P t (Q). (3.2) 

Eq.(3.2) will give very slowly changing or almost flat distribution if we use appropriate 
trial distribution Pt(Q)- We perform Monte Carlo simulation for U(l) gauge theory 
at (3 = 0.0 (strong coupling limit), 1.0, 2.0, 3.0, 4.0 and 5.0 (weak coupling regions), at 
volume V = 100 and 400. 



Table I 



p 


c{p) 




measured C{j3) 
y=100 F=400 
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oo 


5.99 


5.53 


0.5 


8.410 


16.79 






1.0 


12.30 


12.30 


12.4 


12.0 


1.2 


14.38 


11.98 






2.0 


25.82 


12.91 


24.4 


25.8 


2.4 


33.13 


13.80 






3.0 


45.20 


15.07 


40.5 


45.1 


4.0 


66.19 


16.55 


57.8 


66.3 


5.0 


86.87 


17.37 




85.7 


10.0 


186.8 


18.68 






15.0 


285.8 


19.05 






20.0 


384.6 


19.23 






100.0 


1960 


19.60 







Measurements are made 10 4 times for each set. Measured range of Q are from to 
<?Max = 15 and from to 30 for V = 100 (=10 x 10) and 400 (=20 x 20), respectively. 
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P(Q) changes by O(10" 56 ) for j3 = 4.0 at V = 100 and by O(10" 82 ) for (3 = 5.0 at F=400 
between Q = and Q = QMax- 

The distributions obtained through MC simulation show Gaussian distribution of Q 
for all the cases. They are shown in fig. 3 (a) and (b). 



Fig. 


3(a) 




Fig. 


3(b) 



The data for P(Q) is taken by the histogram method and the multinomial distribution 
is assumed to obtain statistical error in each set. The obtained error is ~ 2% of the 
measured P{Q) and much smaller than the symbol of data points in the figure. 

The slope Kv(/3) of Q 2 in the Gaussian form coincides quite well with the analytical 
prediction of eq (2.23). The slope parameter Kv{(3) obtained by numerical calculations 
shows 1/V law for each (3 and the value C(f3) = Kv{(3)V ( eq.(2.23) and (2.24)) is almost 
l/-independent( Table (I)). The fact that all the P(Q)'s measured for /3= 1.0, 2.0, 3.0, 4.0 
and 5.0 have the Gaussian distribution shows that the system undergoes the first order 
phase transition at 6 = tt for all these (3 values. 



4. Conclusions 

We presented the analytic form of partition function Z[9) for U(l) lattice gauge theory 
with 6 -term in 2 dimensions. Character expansion of RSRG can be safely applied to 
the system with #-term. In character expansion method, only orthonormal property of 
irreducible characters is used and the problem of complex Boltzmann factor is harmless 
in this approach. Topological charge distribution for any real coupling constant (3 is given 
analytically, and is shown to have Gaussian form. The coefficient Kv(f3) in Gaussian form 
is analytically given. This form of P(Q) leads to Z(9) expressed by the elliptic theta 
function #3 in V — > 00 limit. This function is known to have infinite numbers of zeros on 
the real negative axis of £ = e ld . The closest zero to — 1 provides the information about the 
deconfining phase transition at 6 = it. The distance from Q = —1 is Im(9 = Ky{f3) oc 1/V 
indicates that the phase transition is consistent with the first order one[20l 



The results for P(Q) given above at any (3 is compared with the numerical simulations. 
In numerical simulations the famous problem of complex Boltzmann factor appears. The 
numerical results thus obtained agree with the analytical results quite well. Topological 
charge distribution P(Q) obtained by numerical simulations obeys Gaussian form. The 
coefficient Kv((3) of Q 2 in the exponent of the Gaussian is proportional to 1/V. We covered 
a wide range in both of Q (up to QMax = 30) and P(Q) (changes by order 10 -82 ). We 
comment here on numerical simulations of P(Q) of CP 1 model |18[]in two dimensions. The 
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results of the model at strong couplings also show Gaussian form, and the coefficient obeys 
1/V law. 

Poisson sum formula leads to an interesting expression eq.(2.30) for the partition 
function. It is given by periodic Gaussian form. Since 1/kv oc V — > large, this expression 
shows phase transition clearly. For mo — 1/2 < a < mo + 1/2 (i.e. mo — n < 9 < mo + n), 
m = mo dominates due to large coefficient in front of (m — a) 2 . For mo + 1/2 < a < 
mo + 3/2(-7T + 2ivmo < < 3tt + 27rmo), m = mo + 1 dominates. From this argument we 
understand the critical value a c = mo + 1/2 (mo = any integer). Poisson sum formula also 
relates $3(1/, 9) to a modular transformed form q &z{i> ,9 )( (2.32 and (2.33)). 
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Table I 



TABLE CAPTIONS 

Coefficients in Gaussian distributions are shown for f3 at 0, 0.5, 1.0, • • • 100. 



FIGURE CAPTIONS 

Fig. 1 ln[C(/3)] vs. ln[/3]. C((3) is approximately linear function of (3 in weak coupling regions. 

Fig. 2 C{j3)/j3 vs. (3. We can see approximately constant behavior at large (3. 

Fig. 3(a) Topological charge distribution log 10 (P(Q)) vs. Q 2 measured by Monte Carlo simu- 
lations for various /3's. The value of Q ranges from to 15 ( =QMax ) and the lattice 
size L x L is 10x10. The inverse coupling constant f3 are taken to be 0.0( square), 1.0( 
diamond), 2.0( circle), 3.0( triangle), 4.0( cross). The lines show fittings according to 
the Gaussian form. 

Fig. 3(b) Topological charge distribution log 10 (P(Q)). Notations are as in Fig. 3(a). <5Max=30 
and L = 20. Data for f3 = 5.0 (nabla) is added. 
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